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Abstract 

Long distance excitation energy transfer between a donor and an acceptor embedded 
in a polymer chain is usually assumed to occur via the Forster mechanism which predicts 
a R~ 6 distance dependence of the transfer rate, where R is the distance between the donor 
and the acceptor. In solution R fluctuates with time. In this work, a Brownian dynamics 
simulation of a polymer chain with Forster energy transfer between the two ends is carried 
out and the time dependence of the survival probability Sp(t) is obtained. The latter can 
be measured by fluorescence resonance energy transfer (FRET) technique which is now 
widely used to study conformations of biopolymers via single molecule spectroscopy. It 
is found that the survival probability is exponential- like when the Forster radius (Rf) 1S 
comparable to the root mean square radius (L) of the polymer chain. The decay is strongly 
non-exponential both for small and large Rf, and also for large kp. Large deviations from 
the Wilemski-Fixman theory is obtained when Rp is significantly different from L. 

1 Introduction 

Fluorescence resonance energy transfer (FRET) is a powerful technique to study many 
aspects of structure and dynamics of polymers and biopolymers. 1-8 In this technique, the 
polymer is doped with a donor and an acceptor at suitable positions along the chain. 
The donor is initially excited optically by a laser light which subsequently transfers its 
energy to the acceptor which is located at a distance R from the donor. In many ap- 
plications, the distance between the donor and the acceptor is fixed, as in the case of 
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rigid biopolymers. 4 ' 5 ' 7,8 However, in many cases of interest, this distance is a fluctuating 
quantity. 1-3 ' 6 For example, the distance between the two ends of a polymer in solution 
executes a Brownian motion. If the polymer is assumed to be an ideal Gaussian chain, 
then the mean square distance between the two ends is Nb 2 , where N is the number 
of the monomer units and b is monomer (or the Kuhn) length. Thus, as the number 
of monomers in the chain increases, the average distance between the two ends also in- 
creases. The probability distribution of R, P(R), however, remains peaked at R = for 
ideal chain, although its height decreases as 1/ y/\N). All these properties play important 
parts in determining the observed dynamics of excitation energy transfer. 

The usually assumed mechanism for excitation energy transfer between a donor and 
an acceptor is the Forster mechanism 6 which gives the following expression for the singlet- 
singlet resonance energy transfer rate, k(R), 

where Rf is the Forster radius and Hf is the rate of excitation transfer when the separation 
between the donor and the acceptor is very small, that is R/Rf — > 0. The Forster radius is 
usually obtained from the overlap of the donor fluorescence with the acceptor absorption 
and several other available parameters 5 . 

The dynamics of Forster energy migration has been traditionally investigated by 
performing time domain measurements of the decay of the fluorescence (due to excita- 
tion transfer) from the donor. 1 ' 4 " 6 More recently, this technique has been used in single 
molecule spectroscopy of biopolymers 7 ' 8 . In the latter, distance dependence of FRET pro- 
vides relevant information about the conformation and dynamics of single biopolymers. 
Recently, FRET from single protein molecules has also been used to study protein folding 9 . 
At any given time after the initial excitation, the fluorescence intensity is a measure of 
the "unreacted" donor concentration. As both /sp and Rp are determined by the donor- 
acceptor pair, the rate of decay of the fluorescence intensity provides a direct probe of the 
conformational dynamics of the polymer. 

When the polymer is in solution, each monomer (or polymer bead) executes Brow- 
nian motion. Because of the connectivity of the polymer chain, this Brownian motion 



of the monomers are strongly correlated. This many body nature of polymer dynamics 
can be described by a joint, time dependent probability distribution P(r N ,t) where 
denotes the positions of all the N polymer beads. The time dependence of the probability 
distribution P(r N ,t) can be described by the following reaction-diffusion equation 10,11 

^P(r N , t) = C B (r N ,t)P(r N , t) - k(R)P(r N ,t) (2) 
where Cp is the full 3N dimensional diffusion operator given by, 

C B (r, t) = Dj2 Ap e?(r , (3 ) 
p[d rj d rj P eq (r,t) 

R is the scalar distance between the two ends of the polymer chain and D is the center of 
mass diffusion coefficient. 

The solution of Eq.2 with the sink term given by the Forster expression is highly non- 
trivial. In two seminal papers, Wilemski and Fixman presented a nearly analytic solution 
of the problem for any arbitrary sink. 10,11 The WF theory has been tested {only for the 
average rate) by computer simulations when the sink is a Heaviside function. 11-13 We are 
not aware of any such simulation study with the Forster rate. Such a study is clearly 
important. 

The objective of this paper is mainly two fold. First, to present results of detailed 
Brownian dynamics (BD) simulations of Eq.2, with k(R) given by the Forster rate (Eq.l). 
Detailed investigation into the time dependence of the survival probability at a time t after 
the initial excitation (Sp(t)) is presented. We believe that this is the first such calculation 
of Sp(t) for this kind of a problem. We find that the time dependence of Sp(t) can be 
non-exponential for a large range of the relevant parameter space (N,/cp, Rf). Second, we 
present a detailed comparison of the simulated rate with the WF theory. This comparison 
has been carried out at the level of time dependent survival probability, again, we believe, 
for the first time. 

The organization of the rest of the paper is as follows. In next section, we present 
the details of the BD simulation. In section III, we discuss the WF theory. In section 
IV we present the simulation results and the comparison with the WF theory. Section V 
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concluded with a discussion. 



2 Simulation Details 

Brownian Dynamics (BD) simulations of polymer motion have been carried out with an 
idealized Rouse chain, in which the set of beads are connected by the harmonic potential, 

o N-l 

^ = ^E fa - r ]+ i? (4) 

where the position vector of a bead j is denoted by rj, and N is the number of beads 
constituting the polymer chain. The mean squared bond length is b 2 . Equilibrium end- 
to-end distance of the polymer chain is, 

((r N -r ) 2 )=L 2 = Nb 2 (5) 

In Rouse model, there exists no excluded volume forces and the hydrodynamic interactions 
between the monomer beads are ignored 13 . More details on this model can be found 
elsewhere 13 . In our present study the polymer chain is additionally characterized by the 
presence of two reactive end groups. This essentially implies that within the time interval 
At, the two end groups react with a probability k(R)At. 

The initial configuration for each trajectory has been selected from a Monte Carlo 
generated random distribution. The following equation of motion has been used in the 
simulations, 

rj(t + At) = rj (t) + Fj(t)At + AX G (t) (6) 

where the positions of the j-th bead at time t and t+At are denoted by rj(t) and rj(t+At), 
respectively. Fj(t) is the total force acting on bead j and AX G (t) is a random Brown- 
ian displacement, taken from a Gaussian distribution with zero mean and the variance 
((X G ) 2 ) = 2At. In writing Eq.l, we have set both Dq and KbT equal to unity; the latter 
is Boltzmann constant times the temperature. Here we adopted a time step, At = 0.001. 
However, smaller time steps as low as 0.0001 has been employed in the limit of both large 
&f and large Rp values, to account for the faster dynamics. During the simulation both 
the mass m and the root mean-square bond length b 2 of the bead have been set to unity 
for computational convenience. 
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The trajectory generated by using the above procedure needs to be terminated when 
the two end groups react. This has been done in simulations in the following way. Each 
time the trajectory is updated, the existing end to end distance R is used in Eq.l to 
calculate the distance dependent rate constant, k(R). We then call a random number 
generator to get a value between zero and unity. If this value is less than k(R)At, then the 
trajectory is terminated. Otherwise, the trajectory is continued. One forms a histogram 
over many such trajectories. This procedure generates an irreversible FRET. 

For each polymer chain constructed randomly, we have equilibrated it for 10,000 time 
steps before switching on the reaction. Subsequently, 50,000 to 1 Lakh trajectories with 
different initial configurations were generated and the survival probability S p (t) was ob- 
tained by averaging over all the trajectories. This procedure was systematically applied 
for the polymer chains containing the beads, N = 20, 50 and 100. 

Before proceeding with the simulations of the Forster transfer, we reproduced the 
results of Pastor, Zwanzig and Szabo (PZS) 12 on the mean first passage time with Heaviside 
sink function of infinite strength. Our simulation results agreed with those of PZS within 
the uncertainity given by PZS. 

3 Wilemski-Fixman Theory 10 11 

Several decades ago Wilemski and Fixman (WF) developed a non-trivial theory for the 
diffusion limited intrachain reaction of a flexible polymer. 10,11 To account for the chemical 
reaction they have added a sink term S to the manybody diffusion equation, just as in 
Eq.2. The WF equation of motion is well-known and we present it below for the sake of 
completion 

^P(r N , t) + £ B P(r N ,t) = -k S(R)P(r N , t). (7) 
In the notation of the present work 

The operator Cb{y N ,t) is given by Eq.3. As already mentioned, the treatment of WF is 
general and can be applied to a reaction with arbitrary sink. 
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Let us define a survival probability S p (t) as the probability that the chain has not 
reacted after time t. S p (t) is then given by, 

S p (t) = J P(r N ,t)dr 1 dr 2 ...dr N (9) 

In order to obtain the survival probability WF made a closer approximation, according to 
which the Laplace transform of Sp(t) can be written as, 

SJ 8 ) = - ^2 (10) 

s s*(l + kD(s)/v eq ) 



where D(s) is defined as, 



roc 

D(s) = / e~ st D(t)dt (11) 
Jo 



which is the Laplace transform of sink-sink time correlation function D(t) defined as, 

D(t) = J d 3 R! j rf 3 R 2 5(Ri)5(R 2 )G(Ri,R 2 ,t) (12) 
The Green function appearing in the above equation is given by, 

g(ri,r 2 ,o=^j { {l _ p2f/2 y*p{ — 2^(1-,*) — ) (13) 

where p(t) is the normalized time correlation function of end-to-end vector (R(O).R(i)) / (R 2 ) 
which can be obtained analytically and is given by the following equation, 

= ^ E ^K-Ajt) (14) 

If we neglect excluded volume and hydrodynamic interactions, A; is given by the following 
expression 10 ' 11 . 

A; = 3D (lTT/Nb) 2 , (15) 

Finally v eq is given as, 

lim D(t) = {v eq f (16) 

Note that u eg is the rate when the distribution of the polymer ends is at equilibrium. Thus, 
v eq gives the initial rate of decay of Sp(t) which will show up as the transient behavior. In 
most cases, the rate of decay should become progressively slower, as the population from 
the sink region decreases as the reaction proceeds. 
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Once the choice of sink function specified, it is straight forward to calculate the survival 
probability by utilizing the above set of equations. WF choice was the heaviside sink 
function. Later, Doi showed that WF method is easy to apply if the heaviside sink function 
is replaced with a Gaussian sink function 14 . Bettizzeti and Perico studied the dependence 
of the rate on the choice of sink function with in the frame work WF theory and supported 
the WF closure approximation 15 . Surprisingly, no analysis of the time dependence of the 
survival probability, Sp(t) has ever been reported. In this study we follow the orignal 
scheme proposed by WF to obtain S p (t) analytically. In doing so we use the Stefest 
algorithm to obtain S p (t) through the Laplace inversion of Eq. 10. 

4 Results and Discussion 

Before discussing the results let us describe the scaling that has been used to compare the 
results obtained by simulation with the theory. In the reduced unit notation adopted in 
simulation, the rate constant has been scaled as /cf = kpb 2 / Dq and the real time has been 
scaled by b 2 /Dq. However, in the original WF theory, time is scaled by 6D/L 2 where D 
is the center of mass diffusion. In the free draining limit, so D = Dq/N and L 2 = Nb 2 . 
Thus, the numbers obtained from WF theory is to be converted to the simulation scaling 
for a comparison of results. The Forster radius is scaled by the bead diameter, b. Another 
important parameter in this problem is the root mean square radius of the polymer as 
this determines the end to end distribution. Although we have carried out simulations for 
N=20, 50 and 100, in this report we shall concentrate mostly on N=50. 

Figures 1 and 2 depict the time dependence of the survival probability S p (t) for two 
different values of Rf, Rf = 1 and Rf = 5, respectively, for a fixed N = 50. In figure 1, Hf 
has been varied from 50 to 0.1 , that is, over two orders of magnitude. The decay remains 
non-exponential over the whole range. In figure 2, fcp has been varied from 10 to 0.1. Here 
the decay is exponential-like. These two figures demonstrate the strong dependence of the 
decay profile of Sp(t) on Rp- Note that the earlier experiments which fitted the quantum 
yield to the Forster expression obtained values which are rather large, comparable to the 
ones shown in figure 2. This could have been due to the use of an equilibrium end-to-end 
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distribution in the fitting, instead of a time dependent probability distribution. In model 
calculations, one usually assumes a small value of Rf ( often in the form of a Heaviside 
sink function). This strong dependence of decay profile on Rp could be potentially useful 
in unravelling mechanism and dynamics of energy transfer. 

It is not difficult to understand the above results qualitatively. For an ideal Gaussian 
chain, the maximum in the probability distribution that the two ends are separated by a 
distance R is located at y/\2N/3)b. For N=50, this value is 5.773b. Therefore, when Rp is 
equal to 5, the decay is facilitated by the presence of a large fraction of the distribution at 
a distance of separation where the transfer rate is large. This can explain the exponential- 
like decay for Rp = 5 (figure 2). However, the situation is completely different for Rp = 1. 
Here the probability of finding a polymer with end to end distance so small is negligible 
and the transfer rate where the bulk of the population is located is very small because of 
the strong R dependence of the Forster transfer rate. Therefore, the decay of the survival 
probability starts slowly (Fig.l) and is determined by the interplay between the diffusion 
and the rate. This explains the shape of figure 1. 

The above discussion also suggests that the shape of the survival probability can 
depend rather strongly on the length of the polymer chain. This is because the Forster 
distance for a given donor-acceptor pair is likely to be independent of the length of the 
polymer chain. But the distribution and also the diffusion rate will be determined by N. 
However, this dependence is not trivial and will be discussed elsewhere. 

In figures 3 and 4, we have compared predictions of the WF theory with the simulations. 
In figure 3, Sp(t) is plotted for two very different values of Uf {Kp = 1 and 10) at Rp = 5 
for N = 50. It is seen that while the agreement is satisfactory at short times for both the 
cases, the same is not true at long times, particularly for the smaller kp. In the latter case 
the simulation also finds a larger non-exponentiality, as discussed later. In figure 4, Sp(t) 
is plotted for for Rp = 1 , kp =1.0 (and N = 50). It is seen that the WF theory breaks 
down in this limit. This is one of the main results of the present study. The agreement 
improves for smaller Kp but becomes worse in the opposite limit. 

We have also simulated both larger and smaller chains. Since the size of the polymer 
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scales with N, it is not possible to compare results for different sizes. In figure 5, we show 
the comparison between the simulation results and the WF theory for N = 100 at Rp = 8 
and kp = 0.1. This is the most favorable parameter space for the WF theory. Although 
the simulated Sp(t) decays somewhat faster, the WF theory prediction is not totally off. 
We have not yet searched for any scaling laws (in N, Rp and L) dependence - work is 
under progress in this direction. 

In figures 6 and 7, we present logarithmic plots of Sp(t) to show the extent of non- 
exponentiality, for different values of Rp and kp. It can be seen that while the decay is 
nearly exponential for small Kp (which is expected), it is strongly non-exponential when 
the value of kp becomes comparable to or larger than the bead diffusion rate, Do/b 2 . 

The inability of the Wilemski-Fixman theory to explain the time dependence of the 
survival probability is surprising. We note that earlier theoretical studies have considered 
only the mean first passage time. In figure 8 we have compared the simulated end-to-end 
vector time correlation function (p(t)) with the slightly approximate expression used by 
WF. The agreement is good, as expected. This agreement improves further for larger N. 
Thus, the failure of WF (as shown in Fig.4) must be due to the closure approximation. 

5 Conclusions 

Use of FRET in single molecule spectroscopy of polymers and biopolymers requires ac- 
curate knowledge of the mechanism of energy transfer, more importantly, the distance 
dependence of the transfer rate. The fluorescence quantum yield can provide only an 
average estimate of the distance between the donor and the acceptor if the mechanism is 
well- understood. This could be sufficient for rigid systems. For many systems of interest, 
for example for understanding the dynamics of protein folding or in the collapse of poly- 
mers, one requires the time dependence of the excitation migration. This will be measured 
in terms of the time dependent survival probability. 

In this work we have presented results of detailed Brownian dynamics simulations of 
Forster energy transfer between the two ends of an ideal Gaussian chain. As noted by 
previous workers 12 , this apparently simple problem is actually highly non-trivial because 
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this is a manybody problem. We have calculated survival probability for a large number of 
values of the transfer rate Hf and the Forster radius, Rf- It is found that while the survival 
probability is exponential- like for small values of kp (compared to the monomer diffusion 
rate, D$/b 2 ) and intermediate Rf, it is strongly non-exponential for small (compared to 
L) R F . 

We have compared the results of the simulation with the well-known theory of Wilem- 
ski and Fixman. It is found that the theory is reliable when the Forster radius Rf is 
comparable to the root mean square radius L of the polymer chain and the transfer rate 
kp is comparable to or smaller than the monomer diffusion rate Dq/o 2 . However, the 
agrement is not at all satisfactory in the limit when Rp is much smaller or larger than L. 

What is the reason for the failure of the WF theory when Rf is substantially different 
from the root mean square radius of the polymer? While it is obvious that the WF closure 
approximation is inadequate in many situations, the exact reason for the failure is not clear. 
In fact, for the nature of the decay curve for large or small Rf can possibly be understood 
even from a one dimensional theory, provided the end to end distance correlation function 
p{t) is given. This problem will be discussed elsewhere, [fill 

Note that the distance dependent rate appears in several other chemical processes, 
like in electron transfer reactions where the rate of transfer is known to show an expo- 
nential distance dependence. It will be interesting to study this problem with the method 
employed here. Another important, long standing problem is the study of reactions in 
realistic polymer chains with excluded volume and hydrodynamic interactions. Work in 
this direction is under progress. 
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Figure captions: 

Figure 1. The survival probability obtained from Brownian dynamics (BD) simulations 
of Eq.2 is plotted against the scaled time for several values of fcp at Rf = 1. The curves 
from top to bottom represent the cases with Uf = 0.1, 1, 10 and 50, respectively. 

Figure 2. The survival probability S p (t), obtained from BD simulations is plotted for 
kF = 0.1, 1 and 10 at Rf = 5. Curves from top to bottom show S p (t) at kp = 0.1, 1 and 
10, respectively. 

Figure 3. BD simulation results have been compared with WF theory at a Forster radius, 
Rf = 5. The upper set shows the case with &p = 1 and the lower set is for kF = 10. In 
both the cases, symbols shows the simulation results while the WF theory predictions are 
represented by the full lines. 

Figure 4. WF theory has been compared with the simulation results at a lower value of 
Rf, namely Rf = 1 and for kp=l. WF theory prediction has been shown by the full line 
while the symbols represent the simulation results. As seen from the figure, WF theory 
seems to break down in this limit. 

Figure 5. The comparison between WF theory and simulation results has been shown 
for a larger polymer chain, N = 100, for Rp = 8 and kF = 0.1. Symbols and the full line 
represents the results of simulation and WF theory, respectively. 

Figure 6. The semilog plot of the survival probability S p (t) which has been obtained 
from simulations, is plotted against the scaled time at kF = 0.1 and Rf = 5. This figure 
shows that the decay is nearly exponential for Rp = 5 and small kp. 

Figure 7. The semilog plot of the survival probability, obtained from simulations is 
plotted against the scaled time at kp = 1.0 and Rp = 1. Highly non-exponential behavior 
of S p (t), is very clear in this limit. 

Figure 8. The end-to-end vector time correlation function pit) is plotted against the 
scaled time for a polymer of mean square length L 2 = 50b 2 . Symbols shows the simulated 
p(t) while the pit) obtained from WF expression is shown by full line. 
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